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Abstract 

The properties of old globular cluster systems in galaxy 
halos are used to infer quantitative constraints on as- 
pects of star formation that are arguably as relevant in 
a present-day context as they were during the protogalac- 
tic epoch. First, the spatial distribution of globulars in 
three large galaxies, together with trends in total clus- 
ter population vs. galaxy luminosity for 97 early-type sys- 
tems plus the halo of the Milky Way, imply that bound 
stellar clusters formed with an essentially universal effi- 
ciency throughout early protogalaxies: by mass, always 
e c i = 0.26% ± 0.05% of star-forming gas was converted 
into globulars rather than halo field stars. That this frac- 
tion is so robust in the face of extreme variations in local 
and global galaxy environment suggests that any parcel 
of gas needs primarily to exceed a relative density thresh- 
old in order to form a bound cluster of stars. Second, it 
is shown that a strict scaling between total binding en- 
ergy, luminosity, and Galactocentric position, Ej, = 7.2 x 
10 39 erg (L/Lq) 2 - 05 (r gc /8kpc) -0 - 4 , is a defining equation 
for a fundamental plane of Galactic globular clusters. The 
characteristics of this plane, which subsumes all other ob- 
servable correlations between the structural parameters 
of globulars, provide a small but complete set of facts 
that must be explained by theories of cluster formation 
and evblution m the Milky Way. It is suggested that the 
Eb(L, r gc ) relation specifically resulted Irom star torma- 
tion efficiencies having been systematically higher inside 
more massive protoglobular gas clumps. 

Key words: galaxies: fundamental parameters - galaxies: 
star clusters - globular clusters: general - stars: formation 



1. Introduction 
The old globular cluster systems (GCSs) in galaxy halos 



large-scale star-formation processes in the early universe. 
But at the same time, most new stars today, whether in 
the Galactic disk or in galaxy mergers and starbursts, are 
born not in isolation but in groups. To be sure, the for- 
mation of a bona fide star cluster is a rare event (only a 
very small fraction of young stellar groups emerge from 
their natal clumps of gas as gravitationally bound units) 
but, as will be argued here, it is one that occurs regularly. 
Thus, individual globular clusters must also be viewed — 
even if in a limit — as the products of a robust mode of 
smaller-scale star formation that has always been viable. 

It should come as no surprise, then, that there is no 
definitive theoretical description of the formation history 
of globular clusters; while some concepts have been identi- 
fied that do seem likely to survive as elements of a correct 
theory in the future, at this point there is simply no model 
that can claim completeness. Detailed discussions of the 
many theories already in the literature may be found in, 
e.g., |Ashman fc Zepf (1998)1 or |Meylan fc Heggie (1997) 



The focus here will instead be on recent progress in ex- 
tracting quantitative and (as nearly as possible) model- 
independent constraints from the data on many GCSs. A 



similarly empirical discussion is given by Harris (2000) 



with an eye mostly to implications for large-scale aspects 



of galaxy formation and evolution (see also Ashman & 
Zepf 1998). In what follows, particular emphasis will be 



placed on applications to the problem of generic star for- 
mation on subgalactic scales. 

A serious concern, when trying to use GCSs in this 
way, is the influence of dynamical evolution on the gross 
properties of cluster systems that have been immersed for 
a Hubble time in the tidal fields of their parent galaxies. 
Two-body relaxation and evaporation, disk- and bulge- 
shock heating, chaotic scattering or disruption by a com- 
pact nucleus, and dynamical friction: all of these pro- 
cesses whittle away at clusters individually and collec- 



tively (Spitzer 1987; Aguilar et al. 1988; Ostrikcr et al 



star fo 'mation in both protogalactic and present-day set- 
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special interest in a conference such as this, as they 1989; 


Capuzzo-Dolcetta 1993 




Murali & Weinberg 1997a; 


potentially to shed light on both local and global Murali & Weinberg 1997b; 


Murali & Weinberg 1997c; 


Ves- 



tings. Because of their great age, their ubiquity in galaxies 
of most any Hubble type, and their impressive spatial ex- 
tent (globulars can be found as far as 100 kpc away from 
the centers of large galaxies) , the integrated properties of 
GCSs clearly had to have been influenced by very general, 



perini 1997). In the case of the Milky Way particularly, the 
net effect is to define a roughly triangular region in mass- 
radius space, within which globulars are predicted to sur- 



vive a 10-Gyr dynamical onslaught (Fall & Rees 1977) — 
the implication being that any clusters born outside such 
a "survival triangle" would have disappeared, taking with 



Proc. of 33rd ESLAB Symp. "Star formation from the small to the large scale" (F. Favata, A. A. Kaas & A. Wilson eds., ESA 
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them vital information on their birth properties. However, 
further investigation shows that Milky Way clusters lo- 
cated at large Galactocentric radii (r gc ^ 3 kpc, roughly 
the effective radius of the bulge) do not fill their expected 



survival triangles (Caputo & Castellani 1984; Gnedin & 



Ostriker 1997). Similarly, evolutionary calculations geared 



to conditions in the giant elliptical M87 ( Murali fc Wcin 
berg 1 )97b| ) suggest that the damage done to an initial 



GCS may be largely confined to within an effective ra- 
dius in that galaxy as well. Thus, there are some features 
of globular clusters, and of GCSs, that are not dictated 
purely by evolution from some completely obscured ini- 
tial conditions, and while care must be taken in identifying 
such observables, the task is not an impossible one. 

Following the suggestion that large fractions of the 
GCSs in ellipticals may have formed in major mergers 
(e.g., Schweizer 1987 ; Ashman fc Zepf 1992 ), and the re- 
lated discovery of young, massive star clusters in systems 



clouds of the Milky Way disk ( [Harris fc Pudritz 1994| ). 
These facts have suggested a general physical picture, de- 
veloped by Harris fc Pudritz (1994)| building o n earlier 
arguments by Larson (1988) and Larson (1993), in which 
handfuls of globulars form in each of many protogalac- 
tic fragments whose gas masses (of order 10 9 M Q ), sizes, 
and internal velocity dispersions correspond to disk GMCs 
scaled up by some three orders of magnitude in mass. Such 
a picture has the obvious appeal of being at least concep- 
tually consistent both with current models of hierarchical 
galaxy formation (there is clearly som e affinity with the 
classic picture of Searle fc Zinn 197£ as well) and with 
the observed pattern of present-day star formation. It also 
has more specific attractions (e.g., the protogalactic frag- 
ments and their protoglobulars are presumed to be sup- 
ported largely by nonthermal mechanisms, thus allowing 
for metal-rich protoglobulars with the correct mass scale) 



but many details remain to be worked out. See McLaugh 



much recent discussion in this field has centered on the in- 
terpretation of (often) bi- and multimodal distributions of 
globular cluster colors (as indicators of integrated metal- 
licity) in relation to their host galaxies' dynamical histo- 
ries (e. g., [Zepf fc Ashman 1993|; [Forbes et al. 1997|; |C6te 
et al. l|998 
hardt 1 



like tt 4 Antennae galaxies ( Whitmore fe Schweizer 1995|), fin fc Pudritz (1996)| (and compare the rather different 



view of Elmcgreen & Efremov 1997) for an attempt at a 
quantitative theory for GCS mass spectra which, with the 



Harris & Pudritz (1994) framework as a backdrop, makes 



99S 



Kisslcr-Patig et al. 1998| ; [Kisslcr-Patig fc Geb 
Cote et al. 200C| ). However, as will be discussed 



further in §p[ mctallicity is completely decoupled from the 
other basic properties of individual globulars in the Milky 
Way; thus, while GCS colors are of interest in the context 
of galaxy formation and chemical evolution, they seem 
to be of marginal importance to the star-formation pro- 
cess itself. In this connection, it is worth noting explicitly 
that literally thousands of genuinely old globulars have 
now been identified, in many galaxies besides our own, 
with colors as red as or redder than solar: high metallic- 
ity has never posed any apparent obstacle to the formation 
of massive star clusters. This cautions against theories of 



explicit use of the emerging links between present-day and 
protogalactic star formation. 

In addition to color and luminosity distributions, pho- 
tometric studies of extragalactic GCSs yield estimates of 
their specific frequencies (total cluster populations, nor- 
malized to the parent galaxy luminosity) and spatial struc- 
tures (number density of globulars as a function of galac- 
tocentric position). These are the focus of where it is 
shown that a large sample of early-type galaxies display 
a common ratio of GCS mass to total mass in stars and 
(hot) gas — an apparently universal efficiency for the for- 
mation of globular clusters (McLaughlin 1999). Finally, ^ 
discusses the binding energies of globular clusters in the 



glob ular formation such a s those of Fall fc Rccs (1985) 



and Murray fc Lin (1993) , which, at least in their current 
form, are able to account only for metal-poor objects. 

Of rather more direct relevance to the small-scale prob- 
lem are the luminosity or mass functions of cluster sys- 
tems. It is well known that the overall range of Galac- 
tic globular cluster masses, m ~ 10 4 -10 6 M , and the 
mean value (m) = 2.4 x 10 5 Mq, are essentially universal 
properties of other GCSs. Traditionally, theories of glob- 
ular cluster formation attempted to explain just this ba- 
sic mass scale (e.g., Peebles & Dicke 1968; Fall & Rees 
1985|)~|More recently, however, attention has turned to the 
full mass spectrum, dN / 'dm (the number of clusters with 
mass m in a single galaxy), which — in the regime usually 
observed, m J> 10 5 M Q — is similar (though not identi- 
cal) from galaxy to galaxy; shows no detectable variations 
with radius in any one system; and is remarkably similar 
to the mass function, dN /dm oc m -17 or so, of the dense 
gas clumps currently forming stars in the giant molecular 



Milky Way (McLaughlin 2000). It is shown that a strong 
correlation between binding energy, total cluster luminos- 
ity, and Galactocentric position is instrumental in defining 
a fundamental plane for Galactic globulars — which, by in- 
corporating every one of their many other structural cor- 
relations, systematically reduces these data to the smallest 
possible set of independent physical constraints for theo- 
ries of cluster formation and evolution. 

2. The Efficiency of Cluster Formation 

The possibility that globular cluster systems can be con- 
nected not only to galaxy formation, but to ongoing star 
formation as well, is suggested by the fact that this lat- 
ter process operates largely in a clustered mode. One dra- 
matic example of this is the situation in the Orion B 
(L1630) molecular cloud, where 96% of a complete sample 
of young stellar objects are physically associated with just 
four d ense clumps of gas, each conta ining > 300 Mq of ma- 
terial ( Lada et al. 1991 ; Lada 1992 ) . More generally, Patej 
|fc Pudritz (1994) note that this is just a result of the dif- 
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ferent power-law slopes in the mass function of molecular 
clumps (as above, dN /dm oc m -17 , so that the largest 
clumps, which always weigh in at 10 2 -10 3 M Q , contain 
most of the star-forming gas mass in any molecular cloud) 
and the stellar initial mass function (dAf/dm oc m~ 2 - 35 , 
putting most of the mass in young stars into ^ 1 Mq ob- 
jects). In more extreme environments, the "super" star 
clusters — luminous, blue, compact associations with inte- 
grated properties roughly consistent with those expected 
of young globulars — found in many merging and starburst 
galaxies can accoun t for as much as ^2 0% of the UV light 
from such systems ( Mcurcr et al. 1995| ). 

Again, however, this is not to say that all, or even 
most, stars are born into true clusters that exist for any 
length of time as systems with negative energy. At some 
point during the collapse and fragmentation of a cluster- 
sized cloud of gas, the massive stars which it has formed 
will expel any remaining gas by the combined action of 
their stellar winds, photoionization, and supernova explo- 
sions. If the cumulative star formation efficiency (SFE) of 
the cloud, r\ = M stars /(M st ar S + M gas ), is below a criti- 
cal threshold when the gas is lost, then the blow-out re- 
moves sufficient energy that the stellar group left behind 
is unbound and disperses into the field. The precise value 
of this threshold depends on details of the dynamics and 
magnetic field in the gas cloud before its self-destruction, 
and on the timescale over which the massive stars dis- 
pel the gas; but various estimates place it in the range 
77 cri t ~ 0.2-0.5 (e g, |Hills 198~0| |Mathieu 1983|; [Elmcgreen 
& Cle mens 1985| ; Verschucrcn 1990| ; and see" Goodwin 1997 
for a discussion of globulars specifically). 

A general theory of star formation must therefore be 
able to anticipate the final cumulative SFE in any sin- 
gle piece of gas with (say) a given mass and density, and 
thereby predict whether or not it will form a bound clus- 
ter. No such theory yet exists. It is possible, however, to 
empirically estimate the probability that a cluster-sized 
cloud of gas is able to achieve r\ > ?7 cr it . This probability — 
or, equivalently, that fraction of an ensemble of massive 
star-forming clouds which manages to produce bound stel- 
lar systems — is referred to here as the efficiency of cluster 
formati on. To get a handle on this for globulars in par- 



ticular, McLaughlin 
fraction 



1999 



works in terms of the 



Ql 



= M init 

— ges 



/M 



init 
gas 



mass 



(1) 



where M™* refers to the total gas supply that was avail- 
able to form stars in a protogalaxy — whether in a mono- 
lithic collapse or a slower assembly of many distinct, sub- 
galactic clumps is unimportant — and Mj^ is the total 
mass of all globulars formed in that gas. The advantage 
of this definition for e c i is that the total mass of a GCS is 
expected to be very well preserved over the course of dy- 
namical evolution in a galaxy; presently observed values 
are reasonable indicators of the initial quantities. This is 
ultimately due to the fact that GCS mass spectra are shal- 



low enough that (again, like the clumps in Galactic molec- 
ular clouds) most of any one system's mass is contained 
in its most massive clusters. But most of the destruction 
processes mentioned in §|l| operate most effectively against 
Zow-mass globulars, which may be lost in great numbers 
(substantially affecting the total GCS population, A/tot) 
while decreasing the integrated M gcs by as little as ~25% 



(McLaughlin 1999) 



Until very recently, it was generally assumed that a 
galaxy's total luminosity, or stellar mass, was an adequate 
stand-in for M™!,*. Thus, the number of globulars per unit 
of halo light has long been taken as a direct tracer of the ef- 
ficiency of globular cluster formation in galaxies. However, 
this approach leads quickly to two interesting problems. 

2.1. GLOBAL AND LOCAL SPECIFIC FREQUENCIES 



Specific frequency was originally defined by Harris & van 
|dcn Bcrgh (1981)| as a global property of galaxies. It is 
nominally the ratio, modulo a convenient normalization, 
of the total GCS population to the total V-band light 
integrated over an entire galaxy: 

S N = Mot x io°- 4 ( M -+ 15 ) = 8.55 x 10 7 (Mat/iv,gai) (2) 

Most subsequent studies of GCSs have therefore estimated 
their total populations and cited S^r-values according to 
equation (|^). It is more useful, however, following the 
discussion just above, to refer to total GCS and stellar 
masses; thus, 



S. 



N 



2500 



(m) 



2.4 x 10 5 M m 



T 



7 Mps L 



-^0 



M s1 



(3) 



for a standard mean globular cluster mass (m) and a rep- 
resentative stellar mass-to-light ratio, Ty iga i, appropriate 
to the cores of large ellipticals (which value is used so as 
not to include any nonbaryonic dark matter in the galaxy 
mass). Note that (m) is not observed to deviate signifi- 
cantly from the Milky Way value, either from galaxy to 
galaxy or from place to place within any one system, but 
that Ty ga j does vary systematically, as a function of lumi- 



nosity, among large ellipticals (e.g., van der Marel 1991). 

Global specific frequencies in a large sample of early- 
type galaxies are shown in Fig. [|. (The GCSs of spirals are 
generally less populous, and often more difficult to iden- 
tify, than those in elliptical systems; thus, the data on late- 
type galaxies are relatively sparse.) The square points cor- 
respond to dwarf ellipticals and spheroidals, some in the 
Local Group (the two faintest objects are the Fornax and 



Sagittarius dwarfs) and others in the Virgo cluster ( Pur 
rell et al. 1996| ; |Miller et al. 1998Q ; filled circles represent 



regular giant ellipticals in a wide range of field and clus- 
ter environments (see, e.g., Harris 1991 and Kissler-Patig 



1997); and open circles stand for the centrally located 



galaxies (which often are also the brightest) i n a large 



number of groups and cluster s ( Blakeslee 1997 ; Blakeslee 
|et al. 19971 ; [Harris et al. 199S| ). 
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Three points are immediately apparent. First, among 
normal gE's an average Sn ~ 5 (roughly, M gcs /M stai s ~ 
0.002) is indicated. Second, the specific frequencies of cen- 
tral galaxies in groups and clusters show a systematic de- 
parture from this "typical" value, increasing strongly to- 
wards brighter galaxy magnitudes. And third, while the 
brightest of the dwarf ellipticals have Sn comparable to 
the giants, the ratio increases towards fainter luminosi- 
ties in these small systems. All in all, Sn ranges over 
more than a factor of 20 in early-type galaxies. If the 
ratio Mg C s/M stars were a good approximation to e c i in 
equation (||) , the implication would be that the basic effi- 
ciency of cluster formation also varied drastically — and in 
a non-monotonic fashion — from galaxy to galaxy. This is 
the first specific frequency problem. Al though it has been 
much discussed in t he literature (see McLaughlin 1999 , 
Elmcgrccn 2000, or Harris 2000 for recent reviews and 
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references), no satisfactory explanation (or prediction) of 
it has ever been advanced. 

Very closely related to this, local specific frequencies 
may be defined at different projected radii within a single 
galaxy, by taking the ratio of its GCS surface (number) 
density profile, N c i(R gc ), and its V-band light intensity, 
iy(-Rgc), normalized as in equation (||). (This is then pro- 
portional to the ratio of cluster and field-star mass den- 
sities, just as in eq. ||.) Beyond an effective radius or so 
(where the effects of dynamical evolution on the GCS are 
presumably minimized), it is found in some galaxies that 
local specific frequencies increa.se outwards; in others, they 
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Figure 2. Second specific frequency problem as seen in 
M87. A distance of 15 Mpc to the galaxy has been as- 
sumed. 



to form at larger distances from the centers of galaxies, in 
gas that was presumably at lower ambient densities and 
pressures. This is the second specific frequency problem. 

Figure H illustrates this situation in M87, the cD galaxy 
at the center of the Virgo Cluster and the first system for 



which the effect was shown convincingly to exist (Har- 



remain l constant. Equivalcntly. some galaxies' CCSs arc 
signifi cantly more extended than their stellar halos. while 



ris 1986). The solid line in the top panel is the galaxy 



light profile, derived from the surface photometry of |dc 



others' ! GCS and field-star distributions trace each other 



accurately on large spatial scales. If the simple assump- 
tion Sn oc e c i were applied here, it would appear to sug- 
gest t hat bound stellar clusters were sometimes (but not 



Vaucouleurs fc Nieto (1978) , and the points trace the pro- 
jected GCS number densities (in units of pc -2 and scaled 
up for a direct comparison with the stellar densities) from 



always [ and for reasons completely unknown) more likely 



the combined data of McLaughlin et al. (1993), McLaugh- 
lin (1995)| , and |Harris (1986)| (see |McLaughlin 1999[ ). It is 
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clear that the radial gradient of N c \ is significantly shal 
lower than that of Iy, leading directly to the strongly 
increasing local specific frequency profile in the bottom 
panel. 

The horizontal line in the bottom of Fig. || marks the 
globally averaged S n for M87 as a whole: 14.1 ± 1.6, three 
times hi gher than "normal" for giant ellipticals ( |Harris et 
al. 1998 ). That is, M87 also suffers from the first specific 



frequency problem. This is one strong hint that the two 
Sn problems are really just different aspects of a single 
basic phenomenon. What this might be has become clear 
only with a homogeneous survey of the centra l galaxies 
(simply BCGs hereafter) in 21 Abcl l clusters by Blakcslce 



1997)| , piakeslee et al. (1997)| , and |Blakeslee (1999) 
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2.2. X-RAY GAS AND A UNIVERSAL e cJ 



Figure 1. First specific frequency problem in 97 early- type 
galaxies. 



Blakeslee has found (see also West et al. 1995) that the 
global specific frequencies of BCGs increase systematically 
with the soft X-ray luminosity (averaged over ^500-kpc 



■5 



scales) of the hot gas in their parent clusters. He then uses 
the details of this correlation to argue that the number of 
globulars in the cores of galaxy clusters scales in direct 
proportion to the total mass there, with one found for ev- 
ery 1-2 x 10 9 Mq of stars, gas, and dark matter. This im- 
plies that the first Sn problem stems (in BCGs, at least) 
from a tendency for brighter galaxies to be "underlumi- 
nous," for the amount of gas and dark matter associated 
with them, rather than overabundant in globular clusters. 
Strictly from a star-formation point of view, however, only 
the baryons are of interest; thus, Harris et al. (1998) sug- 
gest, in essence, that if the dark matter were left out, the 
global mass ratio 

Col = AW(M gas + Mstars) (4) 

might itself be constant, not only among BCGs but in 
other galaxies as wellj^] The first specific frequency prob- 
lem would still arise more or less as Blakeslee suggested, 
with the (observed) larger gas fractions M gas /M st ars in 
brighter galaxies resulting also in a higher global Sn °c 



Mgcs/Mstars- In addition, McLaughlin (1999) notes that 



if the local <T c i — defined in the obvious way as a ratio 
of densities — were also constant as a function of radius 
in galaxies, then the fact that the X-ray gas in ellipti- 
cals tends to be hotter and more spatially extended than 
the stellar distribution could cause the second specific fre- 
quency problem in gas-rich galaxies that also have a high 
global Sn- 

The underlying idea here is, of course, that the present- 
day sum (Mgas+Mstars) should be a better indicator of the 
"initial" gas mass in large galaxies, and that equation (||) 
may therefore be more accurate than Sn as an estimate of 
the true cluster formation efficiency in equation (|l|). Ad- 
mittedly, e c \ is still only an observable approximation to 
e c i; indeed, in a hierarchical universe it may be difficult 
to say precisely what is meant by M™£ in the first place. 
Particularly in a large galaxy which may have accreted gas 
(and stars and globulars) over extended lengths of time, 
the observable e~ c i must be viewed as a mass-weighted av- 
erage over a complex evolutionary history including a po- 
tentially large number of galaxy interactions and many 
discrete star-formation episodes. However, the current ra- 
tio of Mgcs to (Afg as + M s tars) does have an essentially uni- 
versal value in old GCSs — including those in regular gE's 
and early-type dwarfs as well as BCGs — arguing that it 
may be a good reflection after all of the real e c i- 



As 



Blakeslee (1999) discusses at length, his ratio of globu- 

can both be 



lars per unit total mass and the ratio of equation 
constants in BCGs if the baryon fraction in the cores of galaxy 
clusters is also roughly universal (on the order of 10%). But 
as reasonable as it may seem, this possibility is in general un- 
proven, and it is not obvious a priori that the two efficiencies 
are necessarily equivalent. Also, the constancy of Blakeslee's 
ratio has been demonstrated neither globally for objects other 
than BCGs nor locally as a function of position inside any 
one galaxy. The behavior of <T c i in general can therefore not be 
anticipated from Blakeslee's work. 





10 4 




10 3 




10 2 




10 


N 

1 


1 


O 


10 3 


Ph 





10 s 








10 








1 




10 3 




10 2 




1 




1 




Figure 3. Constant local cluster formation efficiency in 
M87, M49, and NGC 1399 (McLaughlin 1999). Bold lines 
in all panels trace the sum of star and gas surface densi- 
ties. Broken vertical lines mark the stellar effective radius 
of each galaxy. 



Figure || shows this first for the local e c \ in M87 and 
in the bright ellipticals M49 (also in the Virgo Cluster) 
and NGC 1399 (BCG in the Fornax Cluster): The ratio 
of projected GCS mass densities (£ c i = (m) N c \) to the 
sum of stellar and gas mass densities (Ty jga i Iy + E gas , 
with Ty, ga i measured separately for each system) is con- 
stant beyond an effective radius in the galaxies.^ The con- 
struction of the individual density profiles is discussed in 
detail by McLaughlin (lggl))! , where a comparison of the 
deprojected quantities is also made, confirming the basic 
result and giving essentially the same numbers for e~ c i. Ev- 
idently, including the gas does alleviate the second specific 
frequency problem, in just the sense suggested above: Sn 
increases with radius in M87 because of a locally vary- 
ing gas-to-star mass ratio, rather than any change in the 
fundamental e" c i; and in the comparatively gas-poor M49 
and NGC 1399, there is no Sn problem in the first place. 
Moreover, the first Sn problem is similarly removed when 
the X-ray gas is taken into account: Although the total 
specific frequencies of M87, M49, and NGC 1399 are sig- 
nificantly different (at 14.1, 4.7, and 6-7), their local GCS 
mass ratios at large radii are consistent with a single value: 
£ci = 0.0026±0.0005 in the mean. Finally, since this is also 



2 The departure of the GCS densities from the stellar profiles 
at smaller radii in M87 and M49 is undoubtedly significant, but 
it is not clear whether this is due to a real decrease in the true 
cluster formation efficiency there, or to dynamical depletion of 
the initial GCS, or to substantial dissipation in the gas that 
formed the field stars (after the globulars were already in place) 
in the innermost regions of the galaxies. 
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Figure 4- Total GCS populations and galaxy luminosities 
for the same systems plotted in Fig. [|. Bold lines are the 
relations predicted by a constant global cluster formation 
efficiency of e G i = 0.0026. Light, dotted lines represent 
constant specific frequencies of Sn = 15, 5, and 1.5. From 
McLaughlin (1999). 



independent of galactocentric radius, it is the same as the 
global cluster formation efficiency in each case. 

These specific examples clearly are consistent with the 
notion of a universal e c i. Figure ^, which essentially re- 
plots the Sn data in Fig. [j], confirms it in detail on a global 
scale. To understand the bold curves in this Figure, note 
that equation (ji|) and the definition A/tot 
10 5 Mq) can be written as 



M gcs /(2.4 x 



A/; o t = 4.17x 10 6 ? cl 1 



Ms. 



M stars y vio 12 AW ' (5) 

The heavy solid line in Fig. |5] is just this equation, with 
(i) £ci = 0.0026 fixed; (ii) V-band galaxy luminosities con- 
verted to stellar masses according to the relation Ty lga i 
6.3 M G Lq 1 (L v , g 



i/10 n Lq) 3 ( |van dcr Marcl 199l| ); and 



(iii) galaxy-wide gas-to-star mass ratios estimated from 
a combination of fundamental-plane scalings and the X- 
ray-optical luminosity correlation: M gas /Af st ars ~ 0.55 x 



(^gai/lO 11 ^©) 1 - 5 ( pVlcLaughlin 199S| j. The steep upturn 



in Atot at high £y, ga i (or the sharp increase in BCG spe- 
cific frequency) is thus due to the fast-growing dominance 
of gas over stars, in the face of a constant GCS mass frac- 
tion. Towards lower luminosities, global gas masses be- 
come negligible and S n decreases steadily because of the 
systematic decrease in Ty iga i for fundamental-plane ellip- 
ticals. As a result, at the luminosity of the Milky Way 
spheroid (disk excluded), an e c \ of 0.0026 also accounts 
for the number of halo (metal-poor) Galactic globulars 
(open square in Fig. ^|; see McLaughlin 1999). 

For the early-type dwarf galaxies with I/y jga i < 2 x 
10 9 Lq, the gas-to-star mass ratio in equation (||) has a 
different meaning: The energy of supernova explosions in a 



single burst of star formation in one of these small galaxies 
may have sufficed to expel all remaining gas from its dark- 
matter well, and while any such gas would, of course, no 
longer be directly observable, a proper estimate of e c i must 
still account for it. The bold, dashed line in Fig. [| is one 
attempt to do this. It represents equation (||) given (i) 
the relation (1 + A /'° s s 7M stars ) ~ (W 2 x 10 9 Lq) 



-0.4 



from the theory of Dekel & Silk (1986); (ii) a constant 
TV, ga i = 2 M Q Lq 1 for the stellar populations; and (iii) 
once again, a fixed t c \ = 0.0026. 

Although scatter remains (at the level of factors of 
~2) in the observed A/tot at any given Lv, ga l in Fig. |I| it 
is important that this is more or less random; the mean 
trends in GCS population as a function of luminosity — 
the essence of the first specific frequency problem — can be 
simply explained if the efficiency of globular cluster for- 
mation were constant to first order. Indeed, deviations in 
Fig. H may reflect the scatter of individual galaxies about 
either the fundamental plane or the Lx~Lb correlation 
used to derive the bold lines there, rather than any sig- 
nificant variations in e c \. This requires further study on a 
case-by-case basis, as does the situation in spirals other 
than the Milky Way. Similarly, there is some indication 
(e.g., Mac Low & Ferrara 199£) that the simple treatment 
of galactic winds (i.e., the model of Dekcl & Silk 1986) 
used to correct for the gas lost from dwarf galaxies may be 
inadequate. At this point, however, it is more than plau- 
sible to assert that globular clusters formed in dE's and 
dSph's as in larger galaxies, always in the same proportion 
to the total mass of gas that was initially on hand. 



One important consequence of this is that the effi- 
ciency of unclustered star formation in protogalaxies could 
not have been universal. In both the faintest dE's and 
the brightest BCGs, globular clusters apparently formed 
in precisely the numbers expected of them, while anoma- 
lously low fractions of the initial gas mass were converted 
into field stars. In the case of the dwarfs specifically, if even 
just the idea of the feedback correction above is basically 
correct then all the globulars had to have formed by the 
time a galactic wind cleared the remaining gas; but this 
must have happened before normal numbers of field stars 
appeared. Thus, the gas which formed bound star clusters 
had to have collapsed more rapidly than that which pro- 
duced unbound groups and associations. This implies that 
it was only those pieces of gas which locally exceeded some 
critical density that were able to attain the cumulative star 
formation efficiency of ij <; 20%-50% required to form a 
bound stellar cluster. In addition to this, the uniformity 
of £ci argues — applying as it does over large ranges of ra- 
dius inside M87, M49, and NGC 1399, and from dwarfs in 
the field to BCGs in the cores of Abell clusters — that the 
probability of realizing such a high SFE depended very 
weakly, if at all, on local or global protogalactic environ- 
ment. Quantitative theories of cluster formation should 
therefore seek to identify a threshold in relative density, 
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Sp/p, that is always exceeded by ~0.26% of the mass fluc- 
tuations in any large body of star- forming gas. 

The "relative" aspect of such a criterion is crucial; 
the GCS data militate strongly against any model rely- 
ing on parameters that are too sensitive to environment. 
One such example is the scenario of Elmegreen & Efremov 
(1997)0 in which the pressure exerted by a diffuse medium 
surrounding a dense clump of gas must exceed a fixed, ab- 
solute value in order to produce a high local r\ and a bound 
stellar cluster. However, since pressures vary by orders of 
magnitude in going from dE's to BCGs, or from large to 
small radii in any one galaxy, this idea seems to imply 
systematic variations in e c \ that are not observed. 

BCGs present a complex problem in larger-scale galaxy 
formation, but it is worth noting that a feedback argument 
like that applied to dwarfs may also be relevant to central 



cluster galaxies like M87 (cf. Harris et al. 1998; McLaugh- 
lin 199^). That is, globulars likely also appeared quickly, 
and in normal numbers, in the densest of star- forming 
clumps (perhaps embedded in dwarf-sized fragments) in 
these deep potential wells. The gas more slowly forming 
field stars could have been virialized thereafter, or moved 
outwards in slow, partial galactic winds. The unused gas 
in this case would have to remain hot to the present day, 
and more or less in the vicinity of the parent galaxy, in or- 
der to appear as the X-ray emitting gas that makes e" c i so 
constant in Fig. |J; but this requirement is certainly con- 
sistent with the BCGs being at the centers of clusters. In 
addition, the feedback in this scenario would have more ef- 
fectively truncated the star formation in the lower-density 
environs at larger galactocentric radii in these very large 
systems, thus giving rise to the second specific frequency 
problem as well. There are other possibilities for BCGs, 
however. It is conceivable, for instance, that their "excess" 
gas and globulars were both produced elsewhere in galaxy 
clusters (in failed dwarfs?) and fell together onto the cen- 
tral galaxies over a long period of time. These questions 
need to be examined in much more detail. 



Finally, McLaughlin (1999) argues that the current ef- 
ficiency of open cluster formation in the Galactic disk is 
also ^0.2-0.4% by mass. This figure is much more uncer- 
tain than it is in GCSs, and it is essentially an instan- 
taneous variant of the time-averaged quantity measured 
for globular clusters. Nevertheless, it clearly suggests that 
whatever quantitative criterion is ultimately required to 
explain e c \ — 0.26% in GCSs may very well prove to be 
of much wider applicability. (One exception may be the 
formation of massive clusters in mergers and starbursts, 
where it has been suggested that e c \ ~ 1-10% [e.g., Zepf 
et al 



with a cumulative efficiency n high enough to produce a 
bound globular cluster. The impressive regularity of this 
occurrence is clearly important, as has just been discussed, 
and its rarity is significant as well: the small value of 
e c i = 0.26% implies that the local SFE in an average bit 
of protogalactic gas was much lower than ?7 cr it ~ 0.2-0.5 
(a fact which is also true of molecular gas in the Galaxy 
today). However, these results say nothing of how an ex- 
treme n comes about in any individual gas clump. This is 
another open problem in star formation generally. Its solu- 
tion requires both an understanding of local star formation 
laws (dp*/dt as a function of p gas ) and a self-consistent 
treatment of feedback on small (~10-100 pc) scales. 

The whole issue is essentially one of energetics in a 
compact, gravitationally bound association of gas and em- 
bedded young stars: When does the combined energy in- 
jected by all the massive stars equal the binding energy of 
whatever gas remains? This point of equality and the cor- 
responding t] can in principle be identified for any given 
star formation law and a set of initial conditions in the 
original gas. The difficulty lies, of course, in deciphering 
what these are; but once this is done, progress will also 
have been made in understanding the probability of ob- 
taining rj > ?7 cr it, i.e., the overall efficiency of cluster for- 
mation in One way to begin addressing this complex 
set of problems empirically is to compare the final binding 
energies of stellar clusters with the initial energies of their 
gaseous progenitors. This is a straightforward exercise for 
the ensemble of globular clusters in the Milky Way. 

baito (1979a) evaluated the binding energies Ef, for 
about 10 bright Galactic globulars, along with a number 
of dwarf and giant ellipticals. He claimed that Eb cx M 1 - 5 
for gE's and globulars alike, while the dwarfs fell system- 
atically below this relation (a fact which he subsequently 
attributed to th e effects of large-sc ale feedback such as 
discussed above f|3aito 1979b|; but cf. |Bcnder et al. 199~2]| ). 



Twenty years later, the data required for a calculation of 
binding energy are available for many more than ten glob- 
ular clusters, and they are of higher quality than those 
available to Saito. 

McLaughlin (2000) computes the binding energies of 
109 "regular" Galactic globulars and 30 objects with post- 
core-collapse (PCC) morphologies. The main assumption 
is that single-mass, isotropic King (1966) models provide 
a complete description of the clusters' internal structures. 
Within this framework, the fundamental definition Eb = 
— (1/2) Jq' Airr 2 p4>dr (with r t the tidal radius of the clus- 
ter) may be written as 



99E; Schweizer 1999], However, this conclusion is E b = 1.663 x 10 41 erg 



very uncertain and requires more careful investigation.) 



3. Globular Cluster Binding Energies 

The focus to this point has been on the frequency with 
which ~10 5 -10 6 Mq clumps of gas were able to form stars 



pc/ \M Q pc" 



£{c) , (6) 



where is the model scale radius (King 1966) Ty.o is the 



core mass-to-light ratio; jo is the central I/-band luminos- 
ity density; and £{c) is a well defined, nonlinear function of 
c = log ( r t /r ), obtained b y numerically integrating King 
models ( [McLaughlin 2000[ ). 
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Values of r , j Q , and c are given for all Milky Way clus- 
ters in the catalogue of Harris (1996)| . However, a deter- 
mination of Ty ; o requires measurements of velocity dis- 
persions ctq, and these are available for only a third of 



the sample, in the compilation of Pryor & Meylan (1993) 
For these objects, application of the King-model relation 
Tv,o = 9(Tq / (AttG 7q jo) gives the results shown in the top 
panel of Fig. || The regular globulars there (the 39 filled 
circles) share a single, constant core mass-to-light ratio: 
(log Ty;o) = 0.16 ± 0.03 in the mean, and the r.m.s. scat- 
ter about this is less than the 1-cr observational error- 
bar shown for log Ty,o- (The results for 17 PCC clusters, 
plotted as open squares, are almost certainly spurious [see 
McLaughlin 2000 1. They are shown for completeness but 
not included in any quantitative analyses here.) This is 



con sistent with separate wor k by Mandushcv ct al. (1991) 



and Pryor fc Meylan (1993 



It can safely be assumed that this same T^o applies 
to all other (non-PCC) Galactic globulars, so that Eb can 
be computed from equation (^) given only ro, jo, and c, 
i.e., on the basis of cluster phot ometry or sta r-count data 
alone. If this is done for the full Harris (1996)| catalogue, a 
very tight correlation between binding energy, total cluster 
luminosity, and Galactocentric position is found: 



E b = 7.2 x 10 J9 erg (L/i ) 2 ' 05 (r gc /8 kpc)" u - 4 , (7) 

with uncertainties of about ±0.1 in the fitted powers on 
L and r gc . This relation is drawn as the line in the middle 
panel of Fig. |{| The r.m.s. scatter of the regular-cluster 
data (filled and open circles) about it is no larger than 
the typical 1-er observational uncertainty on log Eb. 

So far as current data can tell, the constancy of T^o 
and the scaling of Eb with L and r gc a re essentially per- 
fect. Now, in the context of King (1966)| models, any glob- 
ular cluster is fully defined by the specification of just 
four (nominally) independent physical quantities. Given 
the results just presented, it is natural to choose these to 
be log T^o, log Eb, the total log L, and the concentra- 
tion parameter c = log (r t /r ). (Additional factors such 
as Galactocentric position or cluster metallicity are quite 
separate from the model characterization of a cluster, and 
they are to be viewed as external parameters.) But the 
tight empirical constraints on T^o and Eb mean that 
they are not actually "free" in any real sense; in practice, 
Galactic globulars are only a two-parameter family, with 
all internal properties set by log L and c. Equivalently, 
the clusters are confined to a fundamental plane in the 
larger, four-dimensional space of King models available to 
them in principle. The top plots in Fig. [5] are then just two 
edge-on views of this plane. I ts properties are discussed in 
detail by McLaughlin (2000) . where this physical approach 
to it is also compar ed to the more statistical tack taken by 
Djorgovski (1995) , who first claimed its existence, and to 



the different interpretation suggested by Bcllazzini (1998) 
The bottom panel shows the third plot possible in the 
physical cluster "basis" chosen here: concentration vs. to- 
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Figure 5. The fundamental plane of Galactic globular clus- 
ters (after McLaughlin 2000). Top panels are two edge-on 
views; bottom is nearly the face-on view. All correlations 
between any other combinations of cluster observables fol- 
low directly from these three relations between Ty;oi Eb, 
c, L, and r„ r . 
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tal luminosity. (This is close to, but not quite, a face- 
on view of the fundamental plane; see McLaughlin 2000|.) 



these fundamental properties of the Galactic GCS were set 



largely by the cluster formation process (see also Murray 



Although it is not one-to-one like those in the top pan- 



els, there is clearly a dependence of c on log L (see also 
van den Bergh 1994 or Bcllazzini ct al. 1996): roughly, 
c w —0.55 + 0.4 log L, but the scatter about this line ex- 
ceeds the observational errorbar on c. Neither the slope 
nor the normalization of this rough correlation changes 
with Galactocentric position, i.e., the distribution of glob- 
ulars on the fundamental plane is independent of r gc . 

The mean core mass-to-light ratio is also independent 
of Galactocentric radius, and none of the distributions in 
Fig. depend on cluster metallicity. Moreover, since any 
other property of a King-model cluster is known once val- 
ues for Ty.o, Eb, c, and L are given, it follows that all in- 
terdependences between any globular cluster observables 
(and there are many; see, e.g., Djorgovski & Meylan 1994) 
are perforce equivalent to a combination of (i) a constant 
Ty i0 = 1-45 M© Lq 1 ; (ii) equation (Q) for Ef, as a function 
of L and r gc ; (hi) the rough increase of c with L; and (iv) 



fc Lin 1992|; |Bcllazzim ct al. 1996| ; |Vesperini 1997| ). 

Figure 1 finally compares the globular cluster energies 
to estimates for the initial values in their progenitors. This 



generic King-model definitions. McLaughlin (2000) derives 
a complete set of structural and dynamical correlations to 
confirm this basic point: only the quantitative details of 
Fig. |5| — and their insensitivity to metallicity — need be ex- 
plained in any theory of globular cluster formation and 
evolution in our Galaxy. 

It is then important that the Eb(L 7 r gc ) and c(L) cor- 
relatio ns are stronger among clusters outside the Solar 



is done for the specific model of Harris fc Pudritz (1994) 
in which protoglobular clusters are embedded in larger 
protogalactic fragments and have properties analogous to 
those of the dense clumps inside present-day molecular 
clouds (see §|l|). In particular, the column densities of the 
protoclusters are postulated to be independent of mass but 
decreasing with Galactocentric radius: M/ttR 2 ~ 10 3 M© 
pc~ 2 (r gc /8 kpc) -1 , which follows from the clumps being 
in hydrostatic equilibrium and from their parent clouds 
being themselves surrounded in a diffuse medium virial- 
ized in a "background" isothermal potential well with a 
circular velocity of 220 km s _1 . This relation then implies 
E b = GM 2 /R = 4.8 xlO 42 erg(M/M©) 1 - 5 (r gc /8 kpc)~ - 5 , 
which is drawn as the broken line in Fig. ^ By construc- 
tion, this is precisely the mass- energy relation obeyed to- 
day by the massive clumps in molecular clouds in the So- 
lar neighborhood. Intrigu ingly, it is als o the M-E\> scal- 
ing originally claimed by 5aito (1979a) for giant elliptical 
galaxies and (bright) globular clusters. 

The dependence of E\> on r gc in such protoclusters is 
nearly the same as that actually found for the globulars 
today. (It similarly accounts for the observed increase of 
cluster radii with r gc [ [Harris fc Pudritz 1994 ; cf. Murray 



circle ( tilled circles in the plots) than among those within 
it (open circles). Given the relative weakness of dynami- 
cal evolution at such large r gc , this is one indication that 
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Figure 6. Binding energy vs. mass for globulars (points; 
solid line) and their gaseous progenitors (broken line) in 
the Galaxy. Total cluster luminosities are converted to 
masses by applying the constant mass-to-light ratio in- 
dicated. 



& Lin 1992 1 — a trend which is, in fact, equivalent to the 
behavior of Ef, in eq. 0.) In Fig. |, the two are taken for 
convenience to be identical, so that the comparison be- 
tween model and observed binding energies there is valid 
at any given Galactocentric position. 

The difference in the slopes of the two Eb{M) relations 
is significant: The ratio of the initial energy of a gaseous 
clump to the final Ejj of a stellar cluster is unavoidably a 
function of its cumulative star formation efficiency r\\ but 
Fig. ^ shows that this ratio of energies changes system- 
atically as a function of mass, and thus that 77 varied as 
well. Moreover, the fact that the difference between initial 
and final Ef, is largest at the lowest masses implies that 77 
had to have been lower in lower-mass protoglobular s. The 
details of this behavior must rely on the density and ve- 
locity structure of the initial gas; the timescale over which 
feedback expels unused gas; re-expansion of the stars after 
such gas loss; and other such specifics which are model- 
dependent to some extent. The inference on the qualita- 
tive behavior of 77 as a function of protocluster gas mass 
is, however, robust. 

A more quantitative discussion of Fig. || — including 
its implications for the mass function of GCSs, which will 
differ from the mass functions of gaseous protoclusters if 
rj varied systematically from one to the other — has to be 
deferred (McLaughlin, in progress). But this evidence for 
a variable star formation efficiency in protoclusters is itself 
a new target for theoretical attack, most likely through a 



10 



general calculation of star formation and feedback such as 
that described at the beginning of this Section. As was also 
mentioned there, if such models can be made to agree with 
the data in Fig. ^], they will likely also shed new light on 
the empirical efficiency of cluster formation, and perhaps 
on other generic properties of GCSs as well — and, thence, 
on larger-scale processes in galaxy formation. 
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